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Abstract. We provide a precise description of the set of residual boundary 
conditions generated by the self-similar viscous approximation introduced by 
Dafermos et al. We then apply our results, valid for both conservative and non 
conservative systems, to the analysis of the boundary Riemann problem and 
we show that, under appropriate assumptions, the limits of the self-similar and 
the classical vanishing viscosity approximation coincide. We require neither 
genuinely nonlinearity nor linear degeneracy of the characteristic fields. 



1. Introduction and main results 

We are interested in systems of conservation laws in one space dimension 

dtU + d4FiU)]^0, (1.1) 

where U — U{t, x) attains values in R", t > and a; G M arc the scalar independent 
variables and the flux function F : K" is regular. System is strictly 

hyperbolic if the Jacobian matrix DF{U) has n real, distinct eigenvalues 

Ai(C/) <A2 ([/)<••• <A„(t/) (1.2) 

for every U G W\ We refer to the books by Dafermos [TU] and Scrre for an 
exposition of the current state of the theory of hyperbolic conservation laws. 

In the present paper, we characterize the set of residual boundary conditions 
generated by the self-similar viscous approximation introduced by Dafermos [9j et 
al. and we then apply our results, valid for both conservative and non conservative 
systems, to the analysis of the so-called boundary Riemann problem. 

The theory of systems of conservation laws poses several challenges: in general 
classical solutions of Cauchy problems starting out from smooth initial data develop 
discontinuities and break down in finite time. Thus, one needs to look for distri- 
butional solutions. However, distributional solutions are not unique and therefore, 
various admissibility criteria, often motivated by physical considerations, have been 
introduced in the attempt at singling out a unique solution. We refer the reader 
to Dafermos [TOj for an extended discussion on admissibility criteria. The Cauchy 
problem obtained by coupling (jl.ip with the initial data U{0, x) = Uq{x) has been 
studied extensively. Existence and uniqueness results for global in time, admissi- 
ble distributional solutions have been established via the random choice method of 
Glimm [T3] and the front tracking algorithm, see Bressan et al [5] and Holden and 
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Risebro [15] . These results hold true under the assumption that the total variation 
TotVar Uq of the initial data is sufficiently small. If this condition is violated, then 
the oscillation or total variation of the solution may blow up in finite time. Exis- 
tence and uniqueness results for data with large total variation have been achieved 
by imposing additional conditions on the structure of the flux function F. In the 
present paper, we only require that the fiux F satisfies the strict hyperbolicity 
condition ()1.2|) and hence, we focus on data of sufficiently small total variation. 

In view of applications coming from physics, it would be natural to construct 
solutions of by taking the limit e — > 0+ of the viscous approximation 

dtZ' + [FiZ')] - ed, [B{Z')d,Z'] . (1.3) 

Here, B{Z) is an n x n matrix which depends on the physical model under consider- 
ation. Convergence results have been established in specific cases, see Dafermos [10] 
and the references therein. In particular, Bianchini and Bressan [4] proved conver- 
gence in the case when B = I, while the proof of the convergence in the general 
physical case is still a challenging open problem. 

A key tool in the analysis of the system of conservation laws (jl.l|) is the study 
of self-similar solutions in the form U{t,x) = Q{x/t). Indeed, the aforementioned 
random choice and front-tracking schemes are implemented by constructing approx- 
imate solutions that locally have the structure of a self-similar function. Also, the 
analysis of self-similar solutions provides information on both the local (in space- 
time) and the long-time behavior of the solution of a general Cauchy problem. 
Self-similar solutions and the so-called Riemann problem have been extensively 
studied, see in particular the works by Lax [2^, Liu [221 [23], Tzavaras [21] and 
Bianchini [3] in the context of data of small total variation. 

To investigate the structure of self-similar solutions, Dafermos [3], Kalasnikov [TH] 
and Tupciev |26| independently introduced the self-similar viscous approximation 

dtU' + d^[F{U')\ ^etd^[B{W)d^W]. (1.4) 

Indeed, because of the "t" factor in front of the second-order term, the viscous 
approximation (|1.4p admits self-similar solutions U^{t,x) = Q^[x/t) which satisfy 
the ordinary differential equation 

Convergence results for (|1.5p have been established, under suitable conditions im- 
posed on the matrix B, by various authors, see in particular Tzavaras [271 128) . 
Andrcianov [2], Joseph and LcFloch [1^ and the analysis and the references in 
Dafermos [IHl Section 9.8]. 

In this paper, we consider the initial-boundary value problem obtained by focus- 
ing on the domain {t,x) G [0, -|-oo[x [0, -|-cx)[. The initial-boundary value problem 
poses all the challenges of the Cauchy problem that have been mentioned (break- 
down of classical solutions, non uniqueness of distributional solutions), as well as 
additional difficulties due to the presence of the boundary. For instance, consider 
the problem obtained by coupling the viscous approximation (jl.3p with the Cauchy 
and Dirichlct data 

Z'{0,x) = Uo Z'{t,0) = Ub, (1.6) 

where Uq and Ub are given constant states in M". Since we are interested in small 
total variation solutions, we focus on the case when \Uo — Ub\ is sufficiently small. 
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If the matrix B is singular (which is the case for most of the physicahy relevant 
systems), then problem (|1.3p . (|1.6p may be ill-posed, namely admit no solutions. 
However, one can obtain a well-posed problem, at the price of higher technicali- 
ties, by relying on a more complicated formulation of the boundary condition (see 
Bianchini and Spinolo [5]). To simplify the exposition, in the following we assume 
that problem (|1.3p . (|1.6|) is well-posed, but our considerations can be extended to 
the more general case of a suitable class of singular matrices considered in [5]. 

We consider the family of initial-boundary value problems obtained by cou- 
pling (|1.3p and (|1.6p and we assume that Z'^ converge as £ — >■ 0+ to a limit U 
in a suitable topology. Having data (|1.6|) . one expects that the solution U admits 
the self similar form U{t,x) ~ Q{x/t) for some function Q and, hence, one recov- 
ers in the limit a solution to the so-called boundary Riemann problem. Moreover, 
because of the assumption that \Uo — Ub\ is sufficiently small, one expects that the 
trace U = limj_j.Q+ Q{^) is well defined. It should be mentioned that such results 
have been established under additional assumptions, see for instance the analysis 
in Gisclon [12] and in Ancona and Bianchini [1]. We emphasize that, in general, 
one has 

U= lim U(t,x)^ lim QiOy^Ub 

and, also, the limit U of (|1.3p . (|1.6p varies if the matrix B varies, and we refer to 
the articles of Gisclon and Serre [T^ [T3| for these observations. 

To investigate the relation between the boundary data Ub and the so-called 
residual boundary condition U == lim^_^Q+ Q{^), we focus on the non characteristic 
case, namely we assume that, for every U £ M", all the eigenvalues of the Jacobian 
matrix DF(U) are bounded away from zero: 

\i{U) <■■■< Xk{U) < -c < < c < Xk+i{U) <■■■< Xn{U) (1.7) 

for some integer 1 < fc < n — 1 and for a suitable constant c > 0. Then, one expects 
(see again [TJ [12] for rigorous results in specific cases) that there exists a function 
V, the so-called boundary layer, which solves the following system: 

( B{V)V' = F{V)~F{U) 

\ ViO) = Ub lim^^+oo Viy) = U , ^'"^^ 

where V denotes the first derivative of the function V . 

Regarding the self-similar viscous approximation of a boundary Riemann prob- 
lem obtained by coupling equation (jl.4p with the data 

t/"(0,x) = [/o t/"(t,0) = C/b, (1.9) 

Joseph and LeFloch established compactness and convergence results in the case 
when B is the identity matrix in [16j or close to the identity in |17j . Moreover, they 
described the self-similar limit U{t,x) = Q{x/t) and showed that Q has bounded 
total variation. Hence, in particular, the trace lim^_j.o+ Q(^), which we denote by 
U again, is well-defined. Furthermore, they investigated the boundary layer and 
proved in Theorem 4.2 of [16] that there exists a boundary layer V satisfying (jl.Sp . 
Their analysis of the boundary layer involves delicate manipulations of the equations 
that rely on the conservative form of (jl.4p . 

In this paper, we provide a different approach to the analysis of the boundary 
layers of the self-similar viscous approximation (|1.4p , (jl.9p . For simplicity, we focus 
on the case when B = I, but our analysis can be extended to the more general 
case considered in [5] . The main differences of our analysis from the analysis in [16] 
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are the following: first, we rely on completely different techniques since here we 
use center-stable manifold tools coming from the area of dynamical systems in the 
spirit of Bianchini and Bressan [l]. Our approach was inspired by the work of 
Dafcrmos jlO] Section 9.8] concerning the limit of the self-similar viscous approxi- 
mation for the Ricmann problem on the whole real line. However, the presence of 
the boundary amounts for additional challenges which arc tackled by employing the 
Center-Stable Manifold Theorem. As a byproduct, our results apply directly to the 
analysis of non conservative systems, namely we handle the limit of the self-similar 
viscous approximation 

dtU' + A{W)d^W = etd^^W (1.10) 

in the case when the n x n matrix A(U) does not necessarily coincide with the 
Jacobian matrix DF(U) of some flux function F : R" — >• IR". On the other side, 
the limit analysis in [16j applies to both the non characteristic and (under genuine 
nonlinearity assumptions) to the boundary characteristic case, while here we restrict 
to the non characteristic case. 

It should be noted that in the non conservative case the distributional solutions 
of the quasilinear hyperbolic system 

dtU + A{U)d^U 

are not defined and we refer to Dal Maso, LeFloch and Murat [TT] for possible 
definitions of weak solutions. Wc also quote Lc Floch and Tzavaras [21] and the 
numerous references therein for the analysis of the Cauchy problems associated to 
the non conservative systems obtained by taking the limit e — > 0+ of (jl.lOp . 

Before stating our main results, we introduce some notations. We consider sys- 
tem 

V" = A{V)V' 

V{0) = Ub lim,_+oo V{y) = U ^'■''> 
and note that systems (|1.8p and (|l.lip coincide in the case when B = I and 
A{V) = DF{V). Next, we write the ordinary differential equation at the first line 
of (|l.lip as a first order system by setting 

V = W 



W' = A{V)W ^^-^^^ 

and, by relying on the non characteristic condition (jl.7|) and the Stable Manifold 
Theorem (see e.g. Perko [24l Section 2.7]), we conclude that there is a fc-dimensional 
stable manifold M'^{U) C R" x R" which is invariant under (|1.12p and consists 
of the data (Vb,Wo) such that the solution of (|1.12p satisfying V{0) = Vq and 
W{0) ~ Wo has the following asymptotic behavior: 

lim \V{y)-U\^0, lim \W{y)\e''y^^ = 0, 

where c is the same constant as in (|1.7|) . We then define the projection A4y{U) C 
R" of 7W([/) onto the first component as follows 

MvitJ) = {Fsuch that {V,W) £ M'{U) for some W) . (1.13) 

Now, wc state the main result of the paper. 

Theorem 1.1. Let condition (|1.7p hold for some constant c > and a natural 
number k, 1 < k < n — 1. Given a state Uq £ R", there is a sufficiently small 
constant S > and a constant /3 > such that, for any Ub and U satisfying 
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\Ub — Uo\ < P5 and \Ub — U\ < j35, we have that Ub G Aiy{U) if and only if the 
following two properties hold. 

(1) For every sequence {C/"^} C W such that U'^ ^ U as e — > 0'*', there exists 
a sequence of solutions : [0,(5] M" to the problem 

d^Q' 



Q^O) = Ub 



(1.14) 



where I denotes the n x n identity matrix. 
(2) By setting V^{C) :~ Q^{eC), the family {V^} converges to a limit function 
V'^ : [0, +cxd[— > R", as e 0+, uniformly on compact sets. The function 
V'^ solves the problem 



dC 




V"{0) 


^Ub 




lim V"{C) 


= u 


lim 


dV° 

dc 





(1.15) 



Before we proceed, we make some remarks. First, it is obvious that the solutions 
(3^(^) to the ordinary differential equations (|1.14p yield the self-similar solutions 
U^{t,x) of PTTU)) by writing U''{t,x) = Q^^x/t). Second, by relying on the Stable 
Manifold Theorem, it follows immediately that if system (|1.15p admits a solution, 
then Ub G My{U). However, the opposite direction of the implication in the 
statement of the theorem is not trivial. In other words, the proof is mainly devoted 
to showing that from Ub G My{U), then properties (1) and (2) of the theorem 
follow. Third, Theorem ll.ll can be reformulated by stating that the boundary layers 
generated by the self-similar viscous approximation are the same as the boundary 
layers (jl.lip generated by the classical vanishing viscosity approximation. Below, 
we discuss an application of this fact concerning the boundary Riemann problem 
and the corresponding precise result is stated in Section [5] as Proposition 12.11 

Consider the self-similar vanishing viscosity approximation 



etdl^U' 



dtW + A{U')d,U 

U'{t,Q)^Ub (1.16) 
C/^(0,x-)-[/o 

and the classical vanishing viscosity approximation 

dtZ' + A{Z')d,Z' ^ edl^Z' 

Z^t,0) = Ub (1.17) 
Z%0,x) = Uo 

of a Riemann problem. Under the assumption that \Ub — Uo\ is sufficiently small, 
convergence results for (|1.16p and (|1.17p have been established by Joseph and 
LeFloch [ini [TT] and by Ancona and Bianchini [T] , respectively. 
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By combining the analysis in Dafermos Section 9.8] and Theorem II. !( we 
infer that, under appropriate assumptions (for completeness we give the precise 
result in Proposition 12.11 of Section the limits of (jl.l6p and (|1.17l) coincide. 
An analogous result has been proved for conservative systems in Christoforou and 
Spinolo [8] in both the characteristic and the non characteristic boundary case. 
Here we extend this analysis to non conservative systems but we focus on the non 
characteristic case. The reason why these results are not a priori obvious is that 
in general, for initial-boudary value problems, the limit of a viscous approxima- 
tion (|1.3p . (|1.6p depends on the choice of the viscosity matrix B and hence, on the 
type of viscous approximation. 

The exposition of the paper is organized as follows: in Section [2] we compare the 
limits of the classical vanishing viscosity and self-similar viscous approximation of 
a boundary Riemann problem. More precisely, in Subsection 12. H we present exist- 
ing results concerning the Riemann problem established by Dafermos [TU] for the 
self-similar viscous approximation and by Bianchini and Bressan [4] and Ancona 
and Bianchini [1] for the vanishing viscosity approximation. Using these results 
and applying Theorem II. 1) in Subsection 12.21 we show that, under reasonable as- 
sumptions, the limits of ()1.16p and ()1.17p coincide. The proof of Theorem 11.11 is 
provided in Section [31 

2. On the CLASSICAL AND THE SELF-SIMILAR VISCOUS APPROXIMATION OF A 



The aim of this section is to compare the limits of the self-similar (|1.16p and the 
classical vanishing viscosity (jl.l7p approximation of a boundary Riemann problem. 
We first quote existing results related to the wave fan curves corresponding to 
these two different approximations and then we show in Proposition 12.11 that the 
two limits are the same. 

2.1. Preliminary results: construction of the wave fan curves. In [2^, Lax 

studied self-similar, distributional solutions of the system of conservation laws (jl.ip 
by imposing some technical assumptions on the structure of the flux function F. In 
particular, for every i = 1, . . . ,n, he constructed the i-wave fan curve T^{si, U~^) 
passing through a given state S M" with Si € M being the variable parameteriz- 
ing the curve. The curve T'^{si, U^) takes values in R" as Si varies in a neighborhood 
of and for fixed Si small enough, the Riemann problem 



admits a self-similar solution which is either a rarefaction, or a single contact discon- 
tinuity or a single shock satisfying the Lax admissibility condition. This condition 
was introduced in the same paper j20j . Lax's construction was later extended to 
very general cases in the works by Liu [HI [13] , Tzavaras and Bianchini [3] . 

Also, Bianchini and Bressan [4] and Dafermos (TO] Section 9.8] described the 
i-wave fan curve obtained by taking the limit £ — !■ 0^ in the classical vanishing 
viscosity 
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(2.1) 



dtZ' + A{z')d,z' 




(2.2) 
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and in the self-similar viscous approximation (jl.lOp . respectively. We go over the 
constructions in [4] and Section 9.8] in Subsections 12.1.11 and 12.1.31 respec- 
tively, while in Subsection 12.1.21 we quote the description of the limit of (jl.l7p 
in the boundary case provided by Ancona and Bianchini [1]. We note that the 
constructions in [U [TOl [Tj include the case of non conservative systems. 

2.1.1. Wave fan curves induced by the classical vanishing viscosity approximation. 
Here, we describe the i-wave fan curve T^{si, [/+) as constructed in [H Section 14]. 

Given a matrix A{U) satisfying strict hyperbolicity (|1.2p . we denote by Ri{U), 
. . . , Rn{U) a basis of right eigenvectors. 

We first assume s,; > and then, consider the following fixed problem: 



V^{T) 

^ w.»(t) 

e.(r) 

where we used the notation 

/V) 



i?j(Vj(s), a;i(s),5i(s)) ds 



= f{r)-g{T) 



(2.3) 



Xi(yi{s), cJi(s),^i(s)) ds, 



while g denotes the concave envelope of / on the interval [0, s^]. Also, the so-called 
generalized eigenvectors Ri and generalized eigenvalues Xi are defined in [H Section 
14] by considering the travelling waves V to (|2.2p . namely solutions of the system 

( v 



w 



AiV)-^I 



w 



(2.4) 



0. 



By restricting system (j2.4p on a center manifold about the equilibrum point V 
U'^, W = 0, £, ~ Ai(J7+), one eventually determines the system 



V =R,{V,u;,,0^^ 



(2.5) 



where LOi is a scalar unknown and the generalized eigenvector Ri is a suitable 
function taking values in R" and satisfying Ri{U~^ , 0, Ai(?7+)) = i?i(J7+). 

By relying on the Contraction Map Theorem, one can show that problem (|2.3p 
has a unique solution (Vi,Wi,^i) belonging to a suitable metric space and then set 
T'(si, If^) = Vi{si). The construction in the case when < follows similarly by 
taking g in (|2.3p to be the convex envelope of / instead of the concave one. 

2.1.2. Vanishing viscosity limit of a boundary Riemann problem. In [T|, Ancona 
and Bianchini showed that, if \Uq — Ub\ is sufficiently small and the non charac- 
teristic condition (|1.7p holds, then the limit £ — > 0+ of the vanishing viscosity 
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approximation (|1.17p is determined by imposing the condition 

= Quo (S, sk+i , . . . s„) = 0s [S, T^+i (sfe+i , . . . , r"(s„, f/o) . . . )) ■ (2.6) 

Here, S £ EJ^ and (ps is a map parameterizing the stable manifold of (|1.12p about 
the equilibrium point {U,0), i.e. the manifold My{U) in (|1.13p can be expressed 
as follows: 

MviU) = {(l)s{S, U) for some vector S <E R''}. 
Also, T'(si,.) is the curve of admissible states defined by Bianchini and Bressan 
in [4] whose construction is sketched in Subsection 12.1.11 

It can be shown that (|2.6p uniquely determines the values of S, Sk+i, . . . , s„ and 
that the solution Z{t,x) obtained as the e — J> 0+ limit of the classical vanishing 
viscosity approximation (|1.16l) admits the following representation: 

U if x/t < ^k+i{sk+i) 

Vj{T) if x/t = ^j{T), ioTT e]0, Sj[ or ]sj, 0[, j = fc + 1, . . . , n 
^i(sj) if ^j{0)<x/t<^j+i{sj+i) for j = fc+l,...,n-l 
Uo if xA>Cn(0), 

(2.7) 

where (sk+i, . . . , s„) are determined by using (|2.6p . Vj and £,j are the solutions of 
the fixed point problem (|2.3p and the trace U is given by 

C7 = T^+i(s,+i,...,T"(s„,C/o)...) . 

2.1.3. Wave fan curves induced by the self-similar viscous approximation. Now, we 
present the construction of the wave fan curves induced by the self-similar viscous 
approximation as was constructed by Dafermos in [101 Section 9.8] 

To construct the i-wave fan curve (j)i{si, U^) that emanates from some state , 
we consider the system 

satisfied by the self-similar solutions and we decompose the first derivative along 
the right eigenvectors of A{U): 

^ = j2a,mjiQ'm- (2.9) 

By taking a dual basis of left eigenvectors Li{U), . . . , Ln{U), the component aj 
is given by aj = {Lj,Q^), j ~ l,...,n, where (.,.) denotes the standard scalar 
product. 

Substituting (|2.9p into (|1.5p . we arrive at 



= [>^j{V) - 6]aj- +£ P]hiah au 

^ h,l=l 



where the coefficients jSjia are given by 

MQl = -{DL,{Q')Rh{Q'),Ri{Q')) (2.10) 

and DLj denotes the Jacobian matrix of the field Lj . We now perform a change 
of variables by setting ^ = eC and rescale the components aj by setting ujj := eaj 
for every j = 1, . . . n. The self-similar viscous approximation to (|2.8p is denoted 
by in the new variable C, i.e. V^^(C) = Q'^i^O- However, for the convenience 
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of the reader, we abuse the notation, for now, and we drop the index e from 
since the e dependence is clear in what foUows. In Subsection 13. 2| we wih return 
to the original notation V^. Using decomposition (|2.9p . we rewrite system (|2.8p as 
an autonomous first order system 



v 










n 










= e 




e' 


= 0, 





..,n 



(2.11) 



where ' denotes differentiation with respect to (. By linearizing (|2.1ip about the 
equilibium point 



V = U+, wi =0, , 
we get the linear system 



W„=0, ( = X^iU+), £ = 0, 



v 



for every j = 1, 



^ 



(2.12) 



l)2n-f2 



The center subspace Mi of system (|2.12p consists of the vectors (F, w, ^, e) S 
with cjj = for J 7^ J, thus, Mi has dimension n + 3. By the Center Manifold 
Theorem, there is an (n + 3)-dimensional manifold containing all the solutions 
of (|2.1ip that sojourn in a small enough neighbourhood of the equilibrium point 
([/+, 0, Ai ([/+), 0). Such a center manifold is parameterized by Mi and it is tangent 
to the center subspace at the equilibrium. In general, the center manifold of (|2.1ip 
about (J/"*", 0, Ai(J7"'"), 0) is not unique and therefore, in the following, we fix one 
and we denote it by M.i. We refer the reader to the notes by Bressan in [7] for an 
extended discussion about the Center Manifold Theorem. 

As it is shown in [ini Section 9.8], system (|2.1ip ow M.i is equivalent to 



V 
e' 



= (A«(F, w„e,e)-C)w, 

= e 
-0, 



(2.13) 



where the functions R\ and A' take values in R" and R, respectively, and satisfy 

i?,f([/+, 0,A,(t/+),0) =i?,([/+) and A* ([/+, 0, A,(C/+), 0) = A,(t/+). (2.14) 
It should be noted that system (|2.13p consists of (n + 3) equations. 
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To define the value <j)i[si,U^) attained by the z-wave fan curve of admissible 
states emanating from at s = s; > 0, we consider the fixed point problem 

V^[t) =U+ + j w,(s),C,(s),0)ds 
^^i.r) =/»(r)-V(r) (2.15) 

over the interval [0,5^]. Here, 

Jo 

while denotes the concave envelope of /" on the interval [0, Si]. By relying 
on a fixed point argument, it can be shown that (|2.15p admits a unique solution 
(T4(t), a;,i(T), ^i(r)) over [0,Si]. Then, the i-wave fan curve of admissible states is 
defined by setting (j)i{si,U~^) := Vi{si). 

If Si < 0, the construction is similar, except that the function g in p.lSp is the 
convex envelope of / over the interval [si,0]. Moreover, it can be shown that, by 
construction, the curve (j)i{si, [/+) is tangent to the eigenvector i?i(J7+) at Si — 0. 

2.2. Comparison between the self-similar and the classical vanishing vis- 
cosity approximation. Here, we prove that under appropriate conditions on the 
limit obtained by the self-similar viscous approximation (|1.16p . the solution to the 
boundary Riemann problem obtained via the classical vanishing viscosity approx- 
imation (jl.lTp coincides with the one obtained via the self-similar viscous approx- 
imation (jl.l6p . This is not a priori obvious since it is known that, in general, for 
initial-boundary value problems the limit depends on the type of viscous approxi- 
mation. 

The comparison between the two limits is established using Theorem 11.11 and 
the construction of the wave fan curves described in Subsection 12. II The following 
proposition states the result: 

Proposition 2.1. Let Uq G R" and assume that condition (|1.7p holds for an 
appropriate constant c > and for some natural number fc, 1 < A: < n — 1. Then, 
there exists a sufficiently small constant 6 > such that given Ub with \Uo — Ub\ < S, 
the family of self-similar viscous approximations satisfying (|1.16p converges, up 
to subsequences, to a self-similar limit function U(t,x) = Q{x/t). Assume that the 
following two conditions are both satisfied: 

(a) let U denote the trace 

C7= lim QiO, 
5^0+ 

then U satisfies conditions (1) and (2) in the statement of Theorem \l.l\ 

(b) The function U admits the following representation: 

if x/t < £,k+i{sk+i) 

ifx/t~£_j{T) /orrG]0,Sj[ or]sj,0[,j = k + l,...,n 
ifCjiO) < x/t < forj = fc + 1, . . . ,n - 1 

ifx/t>UO), 

(2.16) 

for suitable small values s^+i, . . . , s„ G R, where the functions Vj and 
are solutions of the fixed point problem (|2.15p . 



U(t,x) 
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Then the limit U{t, x) of the self-similar viscous approximation defined by ()1.16p 
coincides with the unique limit Z(t, x) of the classical vanishing viscosity approxi- 
mation defined by (jl.lTp . 

Proof. The convergence result follows from the analysis in Joseph and LeFloch |16) . 
Indeed, the proof of the compactness result in the first part of [16] can be straightfor- 
wardly extended to the non conservative case. We now focus on the characterization 
of the limit. 

Theorem 11.11 states that condition (a) in the statement of Proposition 12.11 is 
equivalent to assuming that Ub G My{U) and hence by relying on conditions (a) 
and (6), we deduce the relation 

Ub = J^UoiS, Sk+i, ■ ■ ■ s„) = (t>s(^S, 0fc+i(.sfe+i, . . . 0„(s„, C/o) ■ ■ • , (2-17) 

where (ps is as in (j2.6p a function parameterizing the stable manifold and the func- 
tion (pj, j = /c + 1, . . . , n is the j~wave fan curve induced by the self-similar viscous 
approximation and is defined by (|2.15p . We note that Jt/^ is a function from R" 
to R" and we claim that it is locally invertible. Indeed, we recall that the stable 
manifold M.^{U) is tangent at the origin to the stable space V'^{U) and the curve 
(f>j{sj, U'^) is tangent at Sj = to the eigenvector Rj{U^) and therefore, the Jaco- 
bian matrix of J-Ug evaluated at the point S ~ 0, Sfc+i = ■ ■ • = = is the n x n 
matrix whose columns are the eigenvectors Ri{Uo), . • . , i?„(C/o). By relying on the 
Local Invertibility Theorem, we conclude that, if \Uo — Ub\ is sufficiently small, then 
relation (|2.17p uniquely determines the values of S, Sk+i, . . . , s„. 

We can now conclude the proof by combining the following two remarks. First, it 
should be noted that, in general, the center manifold about an equilibrium point is 
not unique and, hence, the functions Ri and in (|2.3p . (|2.5p depend on the choice 
of the center manifold. However, as proven in [U Remark 14.2, page 305], the limit 
Z{t, x) that admits representation (|2.7p does not depend on this choice. Second, 
by construction, the solutions of system (|2.13p satisfying the additional condition 
£ = arc solutions of (j2.4l) lying on a center manifold about the equilibrium point 
{V^ ^ U+,W = 6,^ = Xi{U+)). By comparing ((2T5| and ([23]) and using the fact 
that the fixed point of (|2.3p does not depend on the choice of the center manifold, 
we conclude that the function U in (|2.16p coincides with the unique function Z 
in (|2.7p . The proof is complete. ■ 

Remark 2.1. It should be noted that Joseph and LcFloch in [16l Theorem 4.2] 
showed that the limit U of the self-similar approximation (|1.16p satisfies condition 
(a) of Proposition 12 . 1 1 for the conservative case, A{U) = DF{U). 

3. Proof of Theorem 11.11 
In this section, we provide the proof of Theorem 11.11 

By relying on the proof of the Stable Manifold Theorem (see Perko [24l Section 
2.7]) we have that, if there exists a function F° satisfying (jl.lSp . then Ub G A4y{U). 
Establishing the opposite implication amounts to show that, if Ub G A4y{U), then 
properties (1) and (2) of the theorem hold. To accomplish this, we proceed in four 
steps: 

• in Section l3Tl we restrict system ()2.1ip to a center-stable manifold; 

• in Section 13. 2[ we construct the approximating sequence; 
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• ill Section 
limit; 

• in Section 



we establish the convergence of the sequence and study the 



we eventually conclude the proof of Theorem 11.11 

Note that in the following we fix a state U S K" and we exhibit constants cy and 5 
such that, for every Ub satisfying \U — Ub\ < cv5 and Ub G M.y{U), properties (1) 
and (2) in the statement of Theorem ll.ll are satisfied. The values of cy and 5 a priori 
depend on the state U , however, as a matter of fact, the values of these constant 
only depend on some regularity coefficients of suitable functions constructed in 
Section [3?T] It turns out that the values of these functions are uniformly bounded 
on sufficiently small compact sets and hence, the statement of Thcorcm ll.il is valid 
for every [7 in a small enough neighborhood of a given state Uq. 



3.1. Solutions lying on a center-stable manifold. In this subsection, we ana- 
lyze the solutions lying on a center-stable manifold for the boundary layers. 

By linearizing system ()2.1ip about the equilibrium point {V, lu, e) — {U, 0, 0, 0), 
we obtain 



V 



J = l _ 

uj'j = Xj{U)u!j for every j = 1, . . . , n 



(3.1) 



= 



By the non characteristic condition (|1.7p . we have that Xj{U) ^ for every j = 
1, . . . ,ri and therefore, the center-stable space to p.ip is described by the vectors 
(y^, wi, . . . Wn, ^, e) satisfying V e M" and uju+x = 0, . . . , aj„ = 0. This implies that 
its dimension is n 4- A: -I- 2. In the following, to simplify the exposition we use the 
notation W ~ (wi, . . . ,a;„)*. 

Next, we apply the Center-Stable Manifold Theorem and we refer the reader to 
the book by Katok and Hasselblatt [19] for an extended discussion on this subject. 
Here, we only recall the properties needed in this article. By the Center-Stable 
Manifold Theorem, it follows that in a sufficiently small neighbourhood of the 
equilibrium point ([7,0,0,0) there is a so-called center-stable manifold, which is 
locally invariant for system (|2.1ip and it has dimension n -I- fc -I- 2 as the center- 
stable space. In general, the center-stable manifold of a given equilibrium point 
is not unique. However, any center-stable manifold of the equilibrium (J7, 0,0,0) 
for (|2.1ip is always tangent to the center-stable space at (C7, 0,0,0). Moreover, let 
c be the same constant as in (|1.7p . then any center-stable manifold contains all the 
orbits (T^(C)i W^(C)) C(C)i £(0) that are confined in a sufficiently small neighborhood 
of (C7, 0, 0, 0) and satisfy 



lim 



|y(C) - C/| + |Vl^(C)l + 1^(01 + k(C)l exp(-cC/2) = 



By taking the value of & in the statement of Theorem [TTT] sufficiently small, we now 
fix a center-stable manifold M.'^^ defined in the neighborhood of (C7, 0,0,0) given 

by 

\y-V\<b, \W\<b, \i\<b, \e\<b. (3.2) 
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Note that actually M'^'^ depends on C7, but to simplify notations wc do not ex- 
plicitly indicate this dependence. Also, using the notation Wcs = (wi, . . -iOkY, the 
parameterization of A^^* can be chosen in such a way that 

M"' = {{V,W,^,e) such that W = <i>{V,Wc,,^,e)} (3.3) 

for a suitable function $. In particular, if we restrict system ()2.1ip to Ai^'^, then 
our unknowns are the functions V{(), Wcs{C), C(C)i ^(C)- In the following lemma, 
we derive the ODE satisfied by {V, Wcs,^, e). 

Lemma 3.1. Let Ai'^" be the aforementioned center-stable manifold of ()2.1ip . // 
the constant S is sufficiently small, then any solution of ()2.11|) lying on A4'^^ and 
satisfying bounds p.2p is also a solution of the system 



V = * Wcs 



e' =0 



(3.4) 



Here, the function ^ = \E'(y, Wcs, ^, e) takes values in the space of n x k matrices 
and both the functions A = A(V,Wcs,^,£) and T = T{V,Wcs,^,£) take values in 
the space of k x k matrices. Moreover, A(?7, 0,0,0) is the diagonal matrix with 
eigenvalues Xi{U), . . . Xk{U), while T(t/, 0,0,0) is the identity matrix. 

Proof. First, we observe that the vector (V, W — 0,£,, e) belongs to Ai'^" for any 
V, ^ and £. Hence, 

$(y, Wcs = 0, e) = for all V, ^, e, 

where $ is the function given in p.3p . By relying on the regularity of the map $ 
and using p.3|) . we get the expression 

W = ^{V,Wcs,^,e) = ^cs{V,Wcs,(„£)Wcs (3.5) 

for a suitable regular map $cs taking values in the space of the n x k matrices. 
Since the manifold AA'^'^ is tangent to the center-stable space at the equilibrium 
{V,Wcs,£.,£) = ([/, 0,0,0), then the columns of the matrix $cs(C7, 0, 0, 0) are the 
eigenvectors Ri{U), . . . Rk{U): 

$c.s(t^,0,0,0) = (i?i(^)|...|i?fc(i7)). (3.6) 
It should be noted that the first equation of system ()2.11|) can be written as 

V'^ (^Ri{V)\...\RniV))w 
and thus, by substituting (|3.5p above, we get that 

V = '^iV,Wcs,Le)Wcs, (3.7) 

for a suitable matrix ^I* taking values in the space of the n x k matrices. 
Next, we differentiate p.Sp to get 

W' = <i>'csWcs+<fcsW^,. 
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We also have 

Wcs.i, e)Wcs = BV + CW',, + i'd^^csWcs = 

where the functions B = B{V, WcsA, e) and C = C{V, Wcs,^, e) take values in the 
space of n x n and n x k matrices, respectively and satisfy 

B{V,d,^,e) = Onxn C{V,d,^,e)^Onxk for every 1^, e and e. (3.9) 

The exact expressions of B and C arc not relevant here. 

By substituting ()3.8p into the second equation of (|2.1ip . we obtain 

W = <i>csW^, + B^Wcs + CW',, + £(9e$e.)W^cs = 

= (Diag(Ai, . . . , A„) - e/)$c.M^cs + ^^l^cs, 

where the Diag(Ai, . . . , A„) denotes the n x n diagonal matrix with eigenvalues 
Xi{V), . . . , Xn{V) and / stands for the nxn identity matrix. The function F{V, Wcs,£,, £)Wc, 
takes values in the space of the n x n matrices and comes from the quadratic terms 
in the second equation of (|2.1ip . Hence, by choosing a sufRcienly small constant S 
in (13. 2L we have the bound 



\F{V,Wcs,^,e)\<K\Wcs\ (3.11) 
for some suitable constant K > 0. By rewriting p.lOp . we infer 

(C + $,,)W^^, = (Diag(Ai,. . . , A„) ~a?j<;>csWcs+FWcs-B^Wes-eid^^,,)Wc, 

(3.12) 

and by relying on p.6p and p.9p . we deduce that 

(c + $c.,)(t/, 0,0,0) = (i?i(C7)|...|i?fc ([/)), 

hence all the columns of this matrix are linearly independent at ([7,0,0,0). By 
continuity, the columns of the matrix C + computed at {V^Wcs,^-,e) are lin- 
early independent, provided that {V,Wcs,^,£) satisfies (|3.2p for a sufficiently small 
constant 5. Hence, there exists a function L{V, Wcs,£,, e) taking values in the space 
of the k X n matrices and satisfying 

L{V, W,sA, e) (C + <fcs) {V, W,s, t e) = h ■ 
By multiplying p.l2p on the left by L, we arrive at 

where 

A{V, Wcs,^,e) = iDiag(Ai, . . . A„)$cs + LF - LB^ - eL{d^^cs) 

is a function taking values in the space of fc x fc matrices and satisfying 

A(i7, 0, 0, 0) = Diag(Ai(t7), . . . , Xk{U)) . 

To get the above equality, we use properties p.6p . (13.91) and p. lip . Finally, 
T{V,Wcs,£,,£) = L^cs takes values in the space of fc x fc matrices and by rely- 
ing on p.6p and p.9p . we deduce that T(C7, 0,0,0) is the fc x fc identity matrix. 
This completes the proof of Lemma 13.11 ■ 
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3.2. Construction of the approximating sequence. Here, we construct a so- 
lution to the ODE p.4p derived in the previous subsection. Let {If^} C E" be a 
family of data as in the statement of Theorem 11.11 namely ^ U as e — ;> 0+ . 

To begin with, we fix a positive constant e > and a vector Wcs G . Roughly 
speaking, the goal is to construct a solution to (|3.4p . hence lying on A^'^". such 
that Wcs{0) = Wcs, e(C) = £, lim^^+^o ^(C) = and ^(0) = 0. However, in this 
problem a technical difficulty arises since by imposing ^(0) = 0, system (|3.4[1 yields 
^(C) = eC, whereas the manifold M'^'^ is defined for |<^| < 5 (see formula (|3.2p ). To 
overcome this difficulty, we preserve this bound |^| < (5 by restricting to the interval 
C G [0,(5/e] and replacing the condition lim^_j.+oo V(C) = tj^ with V{5/e) = U'^ . 
Also, to simplify the notation, we use e instead of e and we write and W^g to 
emphasize the dependence of the first two components of the solution to (j3.4p on 
the parameter e. After applying the variation of constants formula to p.4|) . we 
obtain the fixed point problem 

W^,{C) = cxp (Ac - ee/2 h) Wcs+ 

+ j^^ exp (a(C -y)- heei'i + hey'' 12 

(3.13) 

where A = A{U, 0, 0, 0). We recall that A is a diagonal matrix having only negative 
eigenvalues and T(i7, 0, 0, 0) = Ik- The following lemma provides existence and 
uniqueness results for the solution of p.l3|) . 

Lemma 3.2. Let and W^g be as in the statement of Lemma \3.1\ If the con- 
stant 6 in the statement of Lemma \3.1\ is sufficiently small, then we can find con- 
stants cv , and cw which do not depend on e and satisfy the following properties: 
< cv,S, cw < 1 o,nd for every sufficiently small e > and for Wcs S with 

\Wcs\<OS, (3.14) 

we have that system (j3.13p admits a unique solution satisfying 

\V%0 -U\<cv6 and |PT/^^(C)| < cu/ cxp ( - cC/2)(5 , (3.15) 

for every C G [0,(5/£]. 

Proof. The proof relies on a standard fixed point argument based on the Con- 
traction Map Theorem. So here, we provide only a sketch of the proof. 
Fix e and Wcs as above and define the spaces 

{V e C°([0, 6/e] ; R") such that \V'iC) - U\ < cyS for every C £ [0, S/e]} 

and 

{W^s G C%[0, d/e] ; M'^) : |M^,^,(C)| < cw<5cxp ( - cC/2)for every C G [0, 5/e]} . 



A 



(v^y),W!giy),ey, 



A 



T{V^iy),W^M,ey,s -I, 



>W^siy)dy, 
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The constants cy and cw are to be determined. By equipping Xy and Xw with 
the norms 

\\V'\\v:=il\\V'\\coao.5M) and \\W^J\w sup cxp (cC/2)|W^,^,(C)|, 

Ce[0,5/e] 

we obtain two closed metric spaces. In the previous expression, 77 denotes a real 
constant to be determined and c is given by (|1.7p . 
Now, we define the map 

: X^r X X^ ^ C°{[0,S/e];m.") x C"([0, (5/e] ; M'=) 

by setting its components Tf (F^, W^J and (V^^, W^^) equal to the first and the 
second component of the right hand side of p.l3p . respectively. 

By direct computations, one can show that it is possible to choose the constants 
Cw, Cy, 6 and rj in such a way that takes values in the set Xy x Xy^ and is, 
moreover, a strict contraction with respect to both the variables and W^^. We 
recall that 9 is the same constant as in (|3.14p . The values of c\y , cy, and 77 depend 
on S, but are independent of e. By applying the Contraction Map Theorem, we 
conclude that ()3.13p admits a unique solution belonging to Xy x Xyy. 

To give a flavour of the construction, we choose to show that 

miVi',W^si)-T!iV,',W;,2)\\v < \{m^Vi\\v + \\W!,,-W!,^\\w) (3.16) 
for {Vf, j) e X|> X X^, i — 1,2. Indeed, we estimate 

\TliV{,W!,,)-T!iVi,W!,,)\iO 

rS/. 

^{V{iy),W^,M,ey,e)^^{Viiy),W;,^{y),ey,e) W^,,iy) 



< 



dy- 



S/e 



* {Vi {y) , W^,, {y) ,ey,e) W^,, {y) - W^,^ (y) 



dy< 



< 



s/e 



M\WU-WlA[v)dy, 



where L denotes a Lipschitz constant of the function with respect to both the 
variables and W^^^, while M is a uniform bound for 4' in the neighborhood of 
([7, 0, 0, 0) given by namely 

\^{V^,Wl,,i,e)\<M , (3.17) 

for every (V^^ Wl^,i,e) such that \V^- - U\<b, \Wl^\ < d, |^| < S and |£| < S. By 
using the estimate 

I W!,i - 1 (y) < exp ( - c y/2) \ \ W^,, - W^,^ \ \w for < y < 5/e, 

we arrive at 



\n{vi,w!,^)-Ti{v^,w, 



= 2//2)||W^c^.i 
2 



cs2)\\v<-[Lcw5\\Vi-Vnv 
which implies p.l6p provided that 5 and rj are sufficiently small. I 



aWw 
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3.3. The structure of the Umit. We now establish the convergence resuU and 
analyze the limit. 



Lemma 3.3. Let the hypotheses of Lemma \3.2\ be satisfied, let S, cy and cw be 

the same constants as in the statement of Lemma \3.S\ and denote by (V^^W^g) the 
solution of system (j3.13p satisfying bounds (j3.15l) . Then, as e — >■ the family 
{{V^ ,W^^)} converges uniformly on compact subsets of [0,+cxd[. The limit V'^ is 
the unique solution of the ordinary differential equation 

{V")" = A{V"){V'^y (3.18) 

satisfying the conditions 

{V°y{Q) = Wcs and lim V°{C) = U (3.19) 

and the bounds 

|(F°)'(C)| <cn/cxp(-cC/2)(5 \V°{0-U\<cvS for every C €[0,+oo[. 

(3.20) 

Moreover, V'^ satisfies 

\V°iO -U\<cvSeiip{-cC/2) for every C e[0,+ool (3.21) 

for a suitable constant cy > depending only on U, cw , c and on the upper bound 
M given by (|3l7| . 

Proof. Wc proceed in three steps. 
Step 1: First, we establish compactness of the sequence {{V^ ,W^g)}. Let M 
denote a constant such that 

\^{V',W!,,^,e)\, \A{V',W!,,te)\, \T{V' ,W!„te) \ < M (3.22) 

for every IV^ — U\ < 6, \W^J < S, |^| < 6 and |e| < S and fix a compact subset [0, 6] 
of [0, +oo[. By relying on p.lSp . we get that satisfies the bounds 

\\V' ~U\\coaoM}<^vS and \\{V'y\\coao,b]} = W'^W^AcoaoM) < McwS (3.23) 

for every < £ < 5/b. Also, W^^ satisfies 

\mjco(m}<^wS (3.24) 

and 

l|(W"-)1lco([o,b]) = ll(A-^CT)^^,^Jco([o,b]) <A^(l + ^)civ'5- (3.25) 

From (|3^ - (|3^ . wc get that the family {{V^,W^^)} satisfies ah the hypotheses 
of Ascoli-Arzela's Theorem. Hence, for any given sequence e„ 0+, there exists a 
subsequence {en^} and a pair {V°,W^g) such that 

V^'-k V° and Wcs" ^ 

uniformly on [0,b] as £„^ 0+. Wc then consider a sequence {bj} of points 
such that bj — > +oo and apply the previous argument on any interval [0,fej]. By 
relying on a standard diagonalization procedure, wc conclude that for any sequence 
£n 0+ there exists a subsequence and a pair {V^, W^^) such that 

V^^k _> yo and Wcs" W° 
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uniformly on the compact subsets of [0, +oo[ as e„j. ^ 0+. In particular, , Wcs" )} 

converges pointwise on [0, +cx)[ and hence from bounds (|3.15p we deduce 

\V°{C) ~U\<cvS and \W°,{C)\ < cw cxp { ~ cC/2)6 for every C £ [0, +oo[. 

(3.26) 

Step 2: Next, we prove that the limit (F", W^^) is independent of the subsequence 
{Suk}- By applying Lebesgue's Dominated Convergence Theorem to p.l3p . we get 
the identities 

T/0(C) ^(v"{y),W^M,0,0)wliy)dy 

J + oo ^ ' 

< 

W^^iO = cxp (AC) + f cxp (A(C - y)) 
^ Jo 

(3.27) 

By relying on a fixed point argument similar to the one in the proof of Lemma l3.2[ 
we get that the solution (y°, of (|3:27| satisfying bounds ([3:26| IS unique. 
Step 3: Last, we establish bound p.2ip . By relying on the representation for- 
mula p.27p and on estimates p.22[) and p.26|) . we deduce that, for every ( S 
[0,+to[, 

I^^O-C^I < ^^^exp(-cC/2)(5:=5yexp(-cC/2)5. 

Finally, by comparing p.27|) with (|2.1ip and p.4p . we obtain that V'^ is a solution 
of the ODE ((XT^ satisfying conditions ([XTO)) . ■ 

3.4. Conclusion of the proof of Theorem ll.il Here, we complete the proof of 
Theorem O 

We first denote by V'^{U) and M'^{U) the stable space and stable manifold, 
respectively, for system (jl.l2p corresponding to the equilibrium point (U, 0). Then, 
we change coordinates using matrix p.6p . More precisely, every point (V, M^) S 
ViU) can be written in the form W = '^csiU, 0, 0, 0)Wcs for some Wcs £ K''- Next, 
we let 9 and S be as in Lemma [3T2] and for given Wcs S satisfying \ Wcs \ < 06, we 
consider the unique solution {V°, W°^) of ([3:27| satisfying bounds ([3:261) . Moreover, 
by varying Wcs, we define the set 

M = [{V%0),<;>cs{U,0,0,0)Wcs) ■■ W^c. eM^ \Wcs\<0S} CW' xW. 

Now, we claim that the set Ai coincides with the stable manifold A^*(C/). To 
prove this claim, we first observe that system (|1.12p is equivalent to system (|2.1ip 
provided that p. lip is restricted to the invariant space {{V,W,^,e) : ^ = £ = 0}C 
R"xR"xMxR. Hence, {(y,W,^,e) : {V,W) G M'{U), ^ = e = 0} C for any 
given center-stable manifold J\A'^'^ of (|2.11l) corresponding to the equilibrium point 
(f7, 0, 0, 0) and we recall that, by Lemma [XTl solutions lying on M'^" satisfy p.4p . 
Thus, to study the stable manifold M.''{U) of (|1.12p . it suffices to study the stable 
manifold to system p.4p . By recalling the proof of the Stable Manifold Theorem 
(see for example Perko [H pp. 104-108]) and the_fact that (F°, W°J satisfies (|3:26)) 
and (|3?27l) we finally estabhsh that M = M^'iU). By taking Ub = ^°(0) and by 
relying on Lemma [3.3) the proof of Part (2) of Theorem 11.11 is complete. Part (1) 
follows easily from Part (2) by defining Q^(^) = V^{£,/e) for £, € [0,6]. 



MV"iy),W°,iy),0,0 



A 



w^s{y)dy 
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